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FOREWORD 


? 


This report has been prepared to expedite early dissemination of the 
information generated under the contract. The data and conclusions must be 
considered preliminary and subject to change as further progress is made on 
this program. This is a progress report covering the work done during the 
third 12 months of the contract; it is not a final report. The NASA Program 
Manager is Dr. C.C. Chamis. 
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1.0 INTRODUCTION 


This technical program is the work of the Engineering Mechanics and Life 
Management Section of GE Aircraft Engines in response to NASA RFP 3-537260, 
"Coupled Structura 1/The rmal/Electromagnetic (CSTEM) Analysis/Tailoring of 
Graded Composite Structures.” The overall objective of this program is to 
develop and verify analysis and tailoring capability for graded composite 
engine structures taking into account the coupling constraints imposed by 
mechanical, thermal, acoustic, and electromagnetic loadings. 

The first problem that will be attacked is the development of plate and 
shell finite elements capable of accurately simulating the structural/thermal/ 
electromagnetic response of graded composite engine structures. Because of 
the wide diversity of engine structures and the magnitudes of the imposed 
loadings, the analysis of these is very difficult and demanding when they are 
composed of isotropic, homogeneous materials. The added complexity of direc- 
tional properties which can vary significantly through the thickness of the 
structures will challenge the state of the art in finite element analysis. We 
are applying AE's 25 years of experience in developing and using structural 
analysis codes and the exceptional expertise of our University consultants 
toward the successful conclusion of this problem. To assist in this, we are 
drawing heavily on previously funded NASA programs. 

We are drawing on NASA programs NAS3-23698, 3D Inelastic Analysis Methods 
For Hot Section Components, and NAS3-23687, Component Specific Modeling, in 
our development work on the plate and shell elements. In addition to these 
two programs, we will draw on NAS3-22767, Engine Structures Modeling Software 
System (ESMOSS) , and NAS3-23272, Burner Liner Thermal/Structural Load Model 
ing, in Task III when we generate a total CSTEM Analysis System around these 
finite elements. This will guarantee that we are using the latest computer 
software technology and will produce an economical, flexible, easy to use 
system. 

In our development of a CSTEM tailoring system, we will build on NASA 
Program NAS3-22525, Structural Tailoring of Engine Blades (STAEBL) and AE 
program, Automatic Improvement of Design (AID), in addition to the program 
system philosophy of ESMOSS. Because of the large number of significant 
parameters and design constraints, this tailoring system will be invaluable 
in promoting the use of graded composite structures. 

All during this program, we will avail ourselves of the experience and 
advice of our Low Observables Technology group. This will be particularly 
true in the Task V proof-of-concept . Their input will be used to assure the 
relevance of the total program. 

Figure 1 shows our program and major contributions in flowchart form. 
This gives a visual presentation to the synergism that will exist between 
this program and other activities. 




Figure 1. Program Flowchart. 
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Figure 2 depicts an integrated analysis of composite structures 
currently under development in the composite users’ community. The severe 
limitations of such a system are not highlighted because three major steps in 
the process are not shown. Figure 3 adds these steps. The analysis system 
really begins with a definition of geometry. A user then defines a finite 
element model simulating this geometry and the anticipated loading. The 
process then moves to defined Step 3._ One cycle through the process ends with 
the prediction of individual ply average stresses and strains. Now comes a 
significant productivity drain, namely, manual intervention to evaluate these 
stresses and strains against strength and durability limits. Based on this, 
the user must decide to (1) change the finite element model, (2) change the 
composite laminate, (3) both of the above, or (4) stop here. 

Obviously, there is a considerable cost savings to be obtained by select- 
ing Number 4. The CSTEM system will obviate the reasons for selecting Number 
4. This system, shown in Figure 4, begins with the definition of geometry, 
as before, but then proceeds to a definition of master regions which contain 
all of the necessary information about geometry, loading, and material proper- 
ties. Step 3 is a constitutive model which develops the necessary structural, 
thermal, and electromagnetic properties based on a micromechanics approach. 
Furthermore, this constitutive model will contain the logic to generate the 
global finite element model based on the variatio n o f the properties, as 
depicted in Figure 5. Using a nonlinear incremental technique, these global 
models will be solved for their structural, thermal, and electromagnetic 
response. Based on this response the global characteristics will be evalu- 
ated, with convergence criteria and decisions made on remodeling. Once the 
global characteristics meet the accuracy requirements, the local characteris- 
tics are interrogated and decisions made on remodeling because of strength, 
durability, or hereditary effects. Once this cycle has been stabilized, 
optimization will be performed based on design constraint. Our goal in Task 
II is to develop finite elements whose characteristics make this system possi- 
ble. Although the structural properties have been highlighted, the thermal 
and electromagnetic properties have as much or more variation, and less work 
has been done in these areas. 


1.1 EXECUTIVE SUMMARF 

"CSTEM” is the acronym for the computer program being developed under 
the NASA contract, "Coupled Structural/Thermal/Electromagnetic Analysis/ 
Tailoring of Graded Composite Structures." The technical objectives for this 
program are to produce radar signal transparent structures having high struc- 
tural performance and low cost. The multidisciplines involved are all highly 
nonlinear. They include anisotropic, large deformation structural analysis, 
anistropic thermal analysis, anisotropic electromagnetic analysis, acoustics, 
and coupled discipline tailoring. The CSTEM system is a computerized multi- 
discipline simulation specialized to the design problems of radar absorbing 
structures. The enabling technical capabilities are implemented in a special 
3D finite element formulated to simultaneously tailor the geometrical, mate- 
rial, loading and environment complexifies of radar transparent structures for 
cost effective optimum performance. 
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Figure 2. Composite Analysis System. 
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Figure 3. Total Composite Analysis System* 
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In each enabling technical discipline a decoupled stand-alone 3D finite 
element code has been developed. An executive program with controlling itera- 
tive solution techniques performs the nonlinear coupling among the particip- 
ating technical disciplines. A geometry and finite element model generator 
specialized for graded composites has been developed as an intimate part of 
this analysis system. 

The structural analyzer is built around the 8-, 16- , and 20-noded iso- 
parametric finite elements with emphasis on the 20 noded. A graded composite 
constitutive model has been developed for these elements which uses the NASA 
Lewis program ICAN as a subroutine to perform composite micromechanics and 
supply, to the analyzer, the requisite composite properties. The composite 
stiffness gradient controls the finite element definition of a structure with 
two major parameters to vary the number of elements through the thickness and 
the number of numerical quadrature points within an element. A unique set of 
local stiffness characteristics is developed for each numerical integration 
point. Integration of these local characteristics over the volume of the 
element provides total element simulation of composite structures including 
such effects as twist-bend coupling. 

The structural analyzer also performs large deformation analysis using a 
unique incremental updated Lagrangian approach with iterative refinement. 
Testing of this capability against classical large deformation problems has 
shown it to be both more accurate and more economical than available alterna- 
tives. Connected with this technical capability is a deformed position eigen- 
analysis capability. All or selected portions of the nonlinear stiffness 
terms can be incorporated into these eigenanalyses . This capability has been 
checked out against available test data and other computer codes . 

In future years of the NASA contract, these capabilities will be combined 
with an optimizer to perform a totally automated integrated analysis of 
graded composite structures. A final task in the program will be a design 
demonstration of the tailoring capabilities of the CSTEM system. 

In order to reach the thermal analysis goals set for CSTEM, the same 8-, 
16- , and 20-noded isoparametric finite elements are utilized. Four heat 
transfer solution options are available: linear steady state, nonlinear 
steady state, linear transient, and nonlinear transient. To overcome the 
previous economic penalty associated with finite element vis-a’-vis finite 
difference heat transfer, a unique solution technique is employed. This is a 
Newton-Raphson iterative technique with right hand side pseudo- fluxes . The 
code will perform the heat transfer analysis of a thermally anisotropic mat- 
erial considering conduction, convection, and radiation. Table I lists the 
parameters involved. 

Routines for the calculation of absorption of electromagnetic waves have 
been written and checked out. The method followed is based on a data bank of 
absorptivity values for given material types specified at discrete values of 
temperatures, frequency, and polarization angle. The information for a 
specific material or materials is read from the data bank file into arrays 


8 


III ll Ih 1 1*4" ■ IHI i I. i I Mi IlllJl :| | lil ||| || III, til I Mill llillllh l| lilm |;||||;| I 



which are used by the program. The absorptivity is linearly interpolated 
from these discrete values to the local values of temperature, frequency and 
polarization angle. 


Table I. Thermal Analyzer. 
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Transient 

Linear 

Nonlinear 

Linear 

Nonlinear 

• 

Temperature 

T 

T 

T(t) 

T(t) 

• 

Time 



t 

t 

♦ 

Thermal Conductivity 


kij(T) 

kij (t) 

kij(T.t) 

• 

Convection Coefficient 


h(T) 

h(t) 

h(T,t) 

• 

Internal Heat Generation 

Qi 

<i 

Q i (t) 

Qi(t) 

• 

Surface Heat Flux 

Q s 

Qs 

Q s ( t) 

Q S (t) 

• 

Convection Boundary 

Q c 

Qc 

Q c (t) 

Q c(t) 

• 

Specified Nodal 

T 

s 

T 

s 

T s (t) 

T s (t) 


Temperatures 





• 

Heat Capacity 

— 

Cp(T) 

— 

C p (T,t) 

• 

Radiation Emissivity 

— 

e(T) 

— 
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2.0 TECHNICAL PROGRESS 


2 . 1 STRUCTURAL ANALYZER 


The CSTEM goals not only require that the finite element be specialized 
for graded composites, but also that it produce accurate large deformation 
analysis and deformed position eigenanalysis both for lower frequency system 
modes and higher frequency acoustic modes. An updated Lagrange large dis- 
placement analysis has been successfully implemented in the CSTEM code. Veri- 
fication and validation were performed using the 20-noded isoparametrics to 
predict the displacements of a pressure loaded cantilever beam. CSTEM results 
correlate very well with published analytical and numerical results as shown 
in Figure 6. For this example, 20 equally spaced nonfollower pressure load 
increments were applied, resulting in a final tip displacement of 7 . 1 inches. 
Appendix A covers the CSTEM large displacement theory. 

Test cases have also been run for the large displacement problem involv- 
ing layered composites using the same 10-inch long, 1-square inch cross 
section cantilever beam model with transverse pressure loads and five aniso- 
tropic layers. Since no analytical results were available for this problem, 
two different finite element models were used to produce verification and 
validation. One model has a single element through the thickness of the beam 
and five elements along the length while the second model has five elements 
through the thickness and five elements along the length. The models and 
comparison of analytical results are shown in Figures 7, 8, and 9. 

To perform the deformed position eigenanalysis, two eigen extraction 
techniques have been incorporated. These are the Determinant Search and 
Subspace Iteration methods with both lumped mass and consistent mass matrices. 
In the determinant search technique, a Sturm sequence check is used to per- 
form an eigenshift to the neighborhood of the desired eigenvalue. Then a 
Rayleigh quotient is employed to iterate the trial vector until convergence is 
achieved. The subspace iteration technique employs a vector transformation to 
reduce the size of the eigenspace. Table II shows an example which highlights 
the benefit of the eigenshift and Rayleigh quotient. 


2.2 THERMAL ANALYZER 

The CSTEM thermal analyzer utilizes the same 8-, 16- , and 20-noded iso- 
parametric finite elements as the structural analyzer. These elements have 
been modified such that they can perform the heat transfer analysis of graded 
composite structures. These structures have in the past presented a large 
obstacle to heat transfer analysis. This difficulty arose because of their 
anisotropic thermal properties and the large gradients in these properties 
that can exist. To overcome these problems, multiple integration points are 
used through the thickness as in the structural element. 

A second major difficulty has to do with the nature of finite element 
heat transfer analysis. In the past, developers have approached the solution 
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Figure 7, 25 Element Model. 



Figure 8. 5 Element Model. 
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Orthotropic Material 
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Figure 9, 10-Inch Cantilever 



Table 11. Eigenshift and Rayleigh Quotient Example. 


L=8 in. 


W»2 in.-*j 




E = 30 x 10 6 psi, V = 0.3, p = 0.298 x 10 s lb/in 3 


Mode 

Shape 

(Axis-About) 

Eigen Values 

Subspace Iterations 

Eigenvalue Shift and 
Rayleigh Quotient 

1st Bend- 
E 

0.357145302 

0.3 + 0.057520891 

1st Bend^ 

0.362975676 

0.4 - 0.037155197 

1st Tors^ 

6.24068457 

6.0 + 0.24067101 

2nd Bend^ 


8.0 + 0.79694164 

2nd Bend^ 

9.53097761 

10.0 - 0.46902631 

1st Axial^ 


16.0 - 0.98479339 

3rd Bendg 


44 + 0.82830848 

3rd Bend^ 

49.4069934 

49 + 0.40698682 
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of finite element heat transfer in a manner similar to the tangent modulus 
technique for material nonlinearity in structural analysis. This has proved 
to be excessively time consuming in computer resources and noncompetitive 
with finite difference approaches. To overcome this problem, the CSTEM 
analyzer uses a right-hand side pseudo flux technique. This strategy not 
only overcomes the computer time penalty but also lends itself to coupled 
solutions. 

The CSTEM thermal analyzer has four solution options: 

• Steady state, linear 

• Steady state, nonlinear 

• Transient, linear 

• Transient, nonlinear. 

In steady-state analysis, the imposed thermal boundary conditions do not 
change with time. In a linear analysis, the thermal properties are constant, 
in a nonlinear analysis the thermal properties are a function of temperature. 
In a transient analysis, the imposed thermal boundary conditions vary with 
time. For a linear transient analysis, the thermal properties are only a 
function of time; for a nonlinear transient analysis, the thermal properties 
are a function of both time and temperature. 

The CSTEM thermal analyzer can perform the heat transfer analysis of a 
three-dimensional graded composite structure considering conduction, convec- 
tion, and radiation with internal heat generation. Imposed boundary condi- 
tions can be nodal temperatures, surface heat fluxes, fluid temperatures, and 
convection coefficients, radiating body temperatures, view factors, and 
emiss ivities . Appendix B is a draft of the theory and User’s Manual. 


2.3 ELECTROMAGNETIC ANALYZER 


The coupled electromagnetic problem for the aircraft engine will be 
treated in two steps: the first is the electromagnetic wave absorption prob- 
lem and the second is the conversion of the electromagnetic energy into heat. 
The engine structure is assumed to be a dielectric material. For any electro- 
magnetic wave propagating through space and impinging on this structure, an 
amount of the electromagnetic energy will be absorbed. The propagation of an 
electromagnetic wave in a dielectric material obeys the following equation. 


2_ /i_ t L. 

3x ye dx 3xJ 3y 





W 2 |J £ <J> = 
r O 0 


0 


( 1 ) 
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where 


- A component of the magnetic field strength vector 5 or a component 
of the electric field vector E 

w = Wave Frequency 

[so = Permeability of free space 

e = Permittivity of free space 
o 

e d = Permittivity of the dielectric 
the relationship between E and H is V X E = - <J-\ H 
Equation (1) may be rewritten as a minimum potential function: 


K*) * l/2/ v [< (ff) + Ky (§*)* ♦ k' (|f ) 2 - A s $ £ ]dV 


where 


K- K r K * ifc- Zd? ife X2 * w2tJ - c 


o o 


Using finite element notation 
[K] e {<J)} e - MM] e {(|»} e = 0 

whe re 


( 2 ) 


[K] e = J v B 1 • [K'] • BdV 


[M] - / v [H.] 1 • [H] dv 


Summing up equation (2) for the entire structure: 


{IK] - X[M]}{4>} = 0 


(3) 
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The solution to ^equation (3), {<(>}, gives the nodal electromagnetic field 

strength vector E or H. Portions of the energy will be absorbed by the 
dielectric structure and converted to heat. The external energy due to the 
electromagnetic field may be written as: 


3W 


— e _ r r K » M sL rM) + K - M <L (M) + k' M <L (4) 

3<t> J v x 3x 34). k 3x ; y 3y 34>. l 3y ; z 3z 3<t>. l 3z J Jav ^ ; 


where 


4 ) = IN. 4 ^ = [H l , H 2 


.H.] 


<h 

4>2 


♦i 


a* 3H i 

M - 7 — 1 
3x i-1 3x 


3H, 

& ■■ 


3H . 


34>. \3x / 3x 


3H. 

K 1 ! M 


Therefore, equation (4) becomes 


§^ } = (/ v B T [K']Bdv]{4») = [K'] e {4>} C 


(5) 


Absorbed energy will convert to internal heat energy as: 


SW^T) 

3T 


[KjHT} 


( 6 ) 


where 

[K x ] = ; v B T • [k] • BdV 

[K^] = Thermal conductivity matrix 

{T} = Vector of nodal temperatures 
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Equating a percentage of (5) to (6) , we get 


[K t ] {T} = a[K fi ] {<*»} 


( 7 ) 


where 


a = the percentage of energy absorbed. 

Equation (7) thus gives the increase in nodal temperature due to 
absorbed electromagnetic energy. 

Routines for the calculation of absorption of electromagnetic waves have 
been written and checked out. The method followed is based on a data bank of 
absorptivity values for given material types specified at d iscret e values of 
temperature, frequency, and polarization angle. The information for a 
specific material or materials is read from the data bank file into arrays 
which are used by the program. The absorptivity is linearly interpolated 
from these discrete values to the local values of temperature, frequency, and 
polarization angle. 

Calculations are made for one given frequency, orientation, and path of 
an electromagnetic wave at a time. Multiple frequencies and/or different 
paths travelled through the structure are handled by separate calculations 
for the different parameters. 

The orientation of the wave is input in the same manner as skew boundary 
conditions where a wave coordinate system is associated with the electro 
magnetic wave. This wave coordinate system is defined such that the direction 
of propagation is along the positive y axis and polarization is measured from 
the positive x axis as shown in Figure 10. Polarization of the wave is 
defined only from 0 to 180 degrees. 

The element face upon which the wave impinges is input and the angle of 
incidence of the wave is calculated from the dot p rodu ct of the inward normal 
of the layer subsurfaces and the positive y axis of the wave coordinate 
system. The angle of incidence must be between 0 and 90 degrees in absolute 
value. The temperature is calculated using the element shape functions. 

The polarization angle is calculated using the dot product between the 
projection of the wave polarization on the layer midsurface and the material 
principle direction. The polarization angle is calculated as being between 0 
and 90 degrees. 

The absorptivity values in the data bank are essentially percentage 
absorption values. The strength of the electromagnetic wave emerging on the 
other side of the layer is calculated as 

Wout = Win * (1.0 - absorptivity) 
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Figure 10. Wave with 90° Polarization. 


where Win is the magnitude of the normal component of the impinging wave. 
Presently, the wave is assumed to pass through the structure on a straight 
path with absorption occurring on the magnitude of the normal component only. 


2.4 CSTEM COUPLED ANALYSIS 


The CSTEM code is a finite element code which has the capability to 
solve structural, thermal, and electromagnetic absorption problems in such a 
way that the effects of one of these phenomena on another are accounted for. 
Some unique features of the code are a 20-noded isoparametric brick with 
multilayer capability, large displacement analysis capability, finite element 
heat transfer analysis with transient capability, generation of material 
properties from constituent properties, an internal composites analyzer, 
internal mesh generation, and calculation of electromagnetic absorption using 
an easily modified data bank of material absorption properties. 

The main executive routine of CSTEM allows the user to perform separate 
structural, thermal, or electromagnetic multiple load case analyses. In 
addition, a combination of any two or all three types of analyses together 
may be done. This is possible due to the modular form of the code. Each 
module is self contained, passing only the required input geometry and control 
information between the modules as well as returning any results which may be 
required as input for an analysis by a following module. 

The structural module uses a 20-noded isoparametric brick and is similar 
to many other isoparametric finite element codes in several ways. It has the 
capability for centrifugal, acceleration, nodal displacement, nodal force, 
temperature, and pressure loadings. The solution technique used is a multi- 
block column solver which allows solution of very large problems since it can 
work on portions of the set of equations separately. 

The more advanced features of the structural module include its ortho- 
tropic material capability. Material properties can be input relative to the 
material axes and then skewed on an element by element basis to obtain the 
desired orientation of the material with the global coordinate system. The 
structural module can also generate the orthotropic material properties it 
needs for composite materials using the constituent properties making up the 
composite. This is done using an internal adaptive version of the computer 
program INHYD, which accesses a data bank containing the material properties 
of the constituents. The properties are calculated based on the volume ratios 
of the constituents . 

Another advanced feature of the CSTEM structural module is its multiple 
layer capability. This allows the modeling of composite structures with many 
material layers without the necessity of using an element for each layer. 
This stiffness of an element with multiple layers is calculated using integra- 
tion points located on the midplane of each layer within the element. The 
stress and strain are then recovered at these same integration points. There 
is a requirement that the element shape follows the layup of the structure so 
that the layers cut through opposite faces of the element at the same height 
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and not diagonally across the element. This requirement points to the use of 
a mesh generator, which is a part of the CSTEM structural module. 

The CSTEM mesh generator is capable of producing various solids of 
revolution from a minimum of input parameters. The generator can produce 
flat surfaces, cylinders, cones, and general double curved surfaces of up to 
360° rotation. At present these different surface types cannot be generated 
together without some manual intervention. 

Another capability that can be used together with the multiple layer 
capability is a composites analyzer, which is adapted from the ICAN computer 
program. This capability must also be used together with the INHYD generation 
of material properties. The composites analyzer takes the stress/strain 
results from the structural module and integrates them through the thickness 
of the structure at some user specified location. This results in a loading 
which can be used by the composites analyzer to do a microanalysis of the 
composite at that particular location. 

The electromagnetic absorption module is based on a data bank of absorp- 
tion properties for different material types. This data bank contains 
absorption properties for the material at discrete values of temperature, 
frequency, and polarization angle. The absorption of electromagnetic energy 
of a specific frequency and polarization by a given material at a specific 
temperature is calculated by linearly interpolating from the discrete data 
bank values. The orientation of an electromagnetic wave is specified similar 
to a skew material so that a coordinate system is associated with the wave 
propagation. This wave coordinate system is defined such that the direction 
of propagation is along the positive Y axis and polarization is measured from 
the positive X axis. The orientation of the wave coordinate system with the 
global coordinate system is specified using skew transformations. 

The element face upon which the electromagnetic wave is impinging is 
specified by the input as well as the path of elements encountered by the 
wave as it passes through the structure. Absorption calculations are made 
for each material encountered and are carried out using midsurface centroid 
values of temperature and orientation. The impingement angle is calculated 
as a dot product of the wave coordinate system Y axis and the midsurface 
centroid normal. The polarization angle is calculated as the dot product of 
the projection on the layer midsurface of the wave polarization and the 
material orientation. 

Absorption calculations are done for one given frequency, orientation, 
and wavepath at a time. If it is necessary to calculate results for several 
frequencies, orientations, or wavepaths, a separate calculation must be done 
for each combination. 

As mentioned these modules may be run separately through the main 
executive routine by the appropriate specification of input variables. 
However, the main intent of CSTEM is to use the results of one module as 
input to another so that the coupling effects of the various phenomena can be 
deduced. The most general procedure followed is to set up a heat transfer 
input deck which describes a time history of heat transfer loadings as well 
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as any temperature dependency of heat transfer properties. A structural 
input deck is then set up with changes in structural loadings with time. 
Structural load case time points must be coincident with heat transfer load 
case time points, although either can contain t ime points between the 
coincident points to account for changes in loadings in one analysis while 
the other remains constant, or to get printout for some intermediate point. 
Finally, a series of electromagnetic absorption analyses are specified at 
each coincident time point if desired. 

The solution proceeds by reading the first structural load case, which 
includes the geometry, structural material properties, and loadings, as well 
as the global control parameters. If a heat transfer analysis is to be done 
along with this particular structural load case the heat transfer module is 
entered and the heat transfer analysis begins by reading the input it 
requires and performing the type of heat transfer analysis requested. The 
heat transfer analysis continues, stepping along in time until it reaches the 
time point where a structural solution is desired. The resultant nodal 
temperatures are passed back to the structural module and the structural 
solution is performed based on the stiffness and loadings at these 
temperatures. After the displacements, stresses and strains are calculated. 
At this time point the analysis of requested cross sections by the composites 
analyzer is done. If any electromagnetic absorption calculations have been 
requested, they are done next using the previously calculated temperatures 
and deformed configuration. This would be the end of a structural load 
case. It may or may not be the end of a heat transfer load case. In either 
case, the next structural load case is read and the process is repeated. This 
continues until the end of the analysis is reached. 
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APPENDIX A 


CSTEM LARGE DISPLACEMENT THEORY 


Large displacement analysis methods can be divided into two types: 
updated Lagrange and total Lagrange. In updated Lagrange large displacement 
analysis the nodal positions are continuously updated to the current equilib- 
rium position and all stress and strain measures are referred to this udpated 
configuration. In contrast, the total Lagrange approach retains the original 
nodal coordinates throughout the analysis with all stress and strain measures 
referred to this original configuration. The udpated Lagrange method handles 
the effects of previous displacements automatically due to the fact that the 
nodal coordinates include these displacements. In the total Lagrange method 
these effects are included by the use of an additional term in the linear 
stiffness and internal force calculations. 

The basic equations are derived from consideration of equilibrium, 
stating that externally applied forces must be balanced by the stresses gener- 
ated internally. For this balance to be correctly achieved mathematically, 
care must be taken to use the proper stress and strain measures. For updated 
Lagrange these are the Cauchy stress and Almansi strain, whereas for total 
Lagrange these are the 2nd Piola-Kirchof f stress and Green strain. 

Both methods use an iterative incremental type of approach. The equilib- 
rium equations to be solved for the updated Lagrange method are 
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are applied loads at time t + At. 
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In this notation, the left superscript refers to the time when the 
quantity is measured or calculated, the left subscript refers to the configu- 
ration that the quantity is referenced to and the superscript in parentheses 
refers to the iteration number. 


The equivalent equations for the total Lagrange method are 


0 w * : [%]) M (i> - * + 1 %“H <t 


-i) 


where 


« H -[ l H* o H o H dv 


o [<] t hJ T l { s } l W dV 

t+ n {■}-/■ : w T o M 

*^°v 

and 

t + At, . are applied loads at time t + At. 

w . 

For the total Lagrange, method an additional term to account for previous dis- 
placements is included in the linear strain-displacement matrix, so this can 
be written as 



The Almansi strain, the strain referenced to the current configuration, 
is used in the updated Lagrange method. The definition of Almansi strain is 


t 

t 




U. . 

i»J 


+ J U. . 

t J,i 


i v 


t 

t 



where tensor notation is used and the left superscript refers to the time 
when the strain occurs and the left subscript refers to the configuration 
that the strain is referenced to. 

The Green strain, the strain referenced to the original configuration, is 
used in the total Lagrange method and is defined as 


24 


Pft if!' i|: iiininiiiii nun mm || | m 



1 e. . = !* ( fc U. . + 1 U. . + 1 u. . U ) 
o ij \o i,j o j ,i o k,i o k,j/ 


These strains are related to the stresses as follows for the Almansi 
strain and Cauchy stress 

t = x C e 

x mn 1 “P* 1 P<1 

and for the Green strain and 2nd Piola-Kirchof f stress 

t S - t c ^ £ 
o ij o ijrs o rs 

These measures with respect to the current and original configurations 
can be related to each other (or referred to any other configuration) through 
the use of the deformation gradient. The deformation gradient is defined in 
matrix notation as 


o M * (o V t { X } T ) T 


where 


3o* 


oV = 


3o* 


3o, 


and 


1 X ) = { ^2 


and T refers to the transpose. 

Using tensor notation the relations between the different measures are 


mn 


= -a ^ t t- 

°p 0 m.i o 0 a *J 
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t e = °X. t e. °X. 
t mn t i,i o ij t j,n 


t p - _£ ^ 

t mnpq Op o m,i o nj o ijrs o p,r o q,s 
where p is the density and 


o 


P 
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APPENDIX B 

CSTEM THERMAL ANALYZER 
DRAFT 

THEORY AND USER'S MANUAL 
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INTRODUCTION 


In recent years the finite element idealization has become a general 
approach for the stress analysis of complex structural system. In order to 
minimize the preparation of data for a thermal stress problem it is desirable 
that the same finite element model be used in the stress and heat transfer 
analysis. For this reason considerable effort is currently being devoted to 
the development of compatible heat transfer and stress analysis programs. 

In general, the stress and heat transfer analysis of solids are coupled. 
The coupling of structural mechanical loads, heat transfer, and electromag- 
netic are discussed in the main programs. Here, we are only talking about the 
heat-transfer portion of the whole program. In heat-transfer analysis non- 
linearities may be due to temperature dependent material properties and, in 
particular, be caused by nonlinear boundary conditions. 

The heat transfer analysis of graded-composite structural problems ideal- 
ized by finite elements is currently in the developmental and experimental 
stage. The procedures of the treatment of temperature dependent material 
properties, and the nonlinear boundary conditions in graded-composite struc- 
tures are covered in the following text. 

The purpose of this report is to explain development of the techniques 
that permit the practical analysis of complex three dimensional heat transfer 
problems of graded-composite structures. 

The Finite Element Formulation 

The three-dimensional isoparametric solid elements are employed. The 

mesh size can be varied and bodies of arbitrary shape can be considered with- 
out difficulty. Material properties can be varied without difficulty. 
Material properties can be different for each element. In general, mixed 
boundary conditions can be handled directly. The equations which govern the 
response of the discrete system generally involve matrices which are symmetric 
and positive definite. Therefore, effective solution techniques can be 
employed for the solution of both the steady state and transient problems. 

An eight, 16, and 20-node, three-dimensional solid element is shown in 
Figure B-l. The natural coordinates (x,y,z) of the eight corner nodes are 
(±1,±1,±1) and of the 12 remaining nodes are (Oil, ±1), (±1,0, ±1), (±1,±1,0). 
The temperature within the element is defined in terms of the nodal 
temperatures by 

6 (x,y,z,t) = e H i (x,y,z) ^(t). 

To express the interpolation functions in a simple form it is convenient to 
define the following 
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Figure B-l. Eight, 16, and 20 Node 
3D Solid Elements. 


Shape Function: 
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where r, s, t are the local coordinates. 
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Equations for Heat Transfer Analysis 

Consider the three-dimensional heat-transfer conditions as shown in 
Figure B-2. In the analysis of heat transfer conditions, we assume that the 
material obeys Fourier's law of heat conduction: 


, 90 ,90 _ . 90 

a* = ' kx j*’ ay = - k y ay» a z - - kz ^ 


where qx, qy, and qz are the heat flows conducted per unit irea, 0 is the 
temperature of the body, and kx, ky, and kz are the thermal conductivities 
corresponding to the principal axes x, y, and z. 

Considering the heat flow equilibrium in the interior of the body we 
thus obtain 


a_ 

8x 


90. 

(kx 


. 3 36. . 3 

* 3^ (ky 3? * ~z 


90. , 

(kz gj) = - ab 


(i) 


where q is the rate of heat generated per unit volume, 
the body the following conditions must be satisfied: • 


On the surfaces of 


0 = 0e 
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Heat Transfer Problems 


convective 3 
heat transfer 

h<W 



C ♦ C ♦€ ♦ C = r 
*1 S 2 S 3 5 4 


specified temperature 


specified heat flow 


Figure B-2. Three-Dimensional Solution Domain for General 
Heat Conduction. 
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where 0e is the environmental temperature, kn is the body thermal conductiv- 
ity, n denotes the direction of the normal to the surface (outward) and £ 
is the heat flow input to the surface of the body. 


Boundary Conditions 

1. Temperature Conditions - The temperature may be prescribed at 
specific points and surfaces of the body, denoted by in 
Equation (2) 

2. Heat Flow Conditions - The heat flow input may be prescribed 

at specific points and surfaces of the body. The heat-flow 

boundary conditions are specified in Equation (3) 

3. Convection Boundary Conditions - Included in Equation (3) are 
convection boundary conditions , where 

q S = h (6 e - 0 S ) (4) 

4. Radiation Boundary Conditions - Radiation g boundary conditions 

are also specified in Equation (3) with q = k (0 - 0 ) 

r (5) 

where 0 is the temperature of the external radiative source and K is a 

coefficient : 

K = h [0 2 + (0 s) 2 ] [0 + 0 S J (6) 

r r r J 

The variable h is determined from the Stefan-Boltzman constant, the 
emissivity of the radiant and absorbing materials and the geometric view fac- 
tors. In addition to these boundary conditions also the initial conditions 
must be specified in a transient analysis. 


Solution Scheme 

For the development of a finite element scheme either a Galerkin formu- 
lation operating on the differential equation of equilibrium or a variational 
formulation of the heat-transfer problem can be employed. In the variational 
formulation a function II is defined such that when invoking the stationarity 
of II, the governing differential Equation (1) to (3) are obtained. 


33 


II = / 1/2 {kx (||) 2 + ky (|V + kz (||) } dv 

V 

- / 0 q b dv - -^ 2 0 S £ S ds - 

where 0* are the concentrated heat flow inputs. Using the condition of 
stationarity of II we obtain 

/ 60't k 6 dv = / 60 £ b dv + f 60 s £ S ds 

v V SZ 


+ I 6 0 1 Q 1 


( 8 ) 


For a general solution scheme of both linear and nonlinear, steady-state 
and transient problems we aim to develop incremental equilibrium equations. 
The modified Newton-Raphson iteration for heat flow equilibrium is used, in 
which 


t+At 


0 


(i) _ 


t+At 


(i-i) 


+ A 0 


(i) 


With the initial condition 0^ = 0 t and A 0 (l) is the temperature increment. 

In the iteration procedure, the stiffness matrix [K] is constant and ^ny 
change is reflected by an update of the right hand side lo^d vectors [Q ] 
at each iteration. Also the "COLSOL" technique is used to solve the matrix 
for the nodal temperature. 

There are four generic heat transfer problems, they are discussed in the 
following. 


I) Linear Steady- State Conditions 

The governing equation of linear steady-state can be written as: 


(K k + 


K C ) 


t+At e = t+At Q + t+At Q e 


where K 
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is the 


conductivity matrix 
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m 



f B^ T . . dV 


(A) 


(I) 


34 



Q 

K is the convection matrix 


C = J, f h (m) . H 8(m)T . H s(m) . ds 

m-l 

sc 


Th“ nodal point heat flow Input vector t+ t} 
t+At Q = t+it q B ♦ t+at Q s + t+it Q c • 


f jj(m)T t+At^b(m) dv 
v 

_ „ s(m)T t+At s(m) , 

J n . 2 . ds 

S 2 

Where ‘'‘“‘"Q is a vector of concentrated nodal point heat flow 

The nodal point heat flow contribution is due to 

boundary condition e 

gb - the rate of heat generated in element 

g 

<1 - the surface heat flow input 

At that time: 
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where (e) denotes element m 


(la) 

(lb) 

(lc) 


(II) 


(B) 


(III) 


(IV) 


input . 

the convection 


35 



t+At 0 a vector of all nodal point temperatures at t+At 


t+At„ [ t+At t+At„ 

0 - 0j e 2 


t+At 0.1 

1 


fe") 

The matrix H + Element Temperature 

g( e ) + Temperature gradient interpolation matrice 

gS(e) ^ fjjg sur face temperature interpolation matrix 

[B] - derivative of the shape function wi|h respect to 
r, s, t and premultiplication by J 


t+At„e 


Q e = j r , ,h ■ H stm) . H s(m) . t+At e . d S <"> 
x m J (m) e 

sc 


(V) 


0 The given nodal point environmental temperatures . 

e t+At 

from this equation we find the t+ ^ t Q e (i.e* and h are given) 


<j - the rate of heat generated per unit volume 


0 - the environmental temperature 

e 

q s - the heat flow input to the surface of the body 

h - the convection coefficient 

q S = h(0 - 0 s ) 

e 

0 - is the temperature of the body 

a) From Equations (I), (II) we can find the 

(K k + K C ) 

b) From Equation (B) the nodal point heat flow input 


t+4t Q = ttat Q B * t+it Q s ♦ t+4t Q c 
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Q c is a vector of concentrated nodal point heat flow input 


It can be solved for 


c) From Equation (V) the given nodal point environment temperatures are 

t+it e 


t+At 


m 


^(nOT „s(m) t+At c 


Q. = 1 - f h . h . ds 

m-1 ^ 


sc 


We can solve for ., A4prt e / , - . ^ x 

t+AtQ (nodal point heat distribution). 

From Items (a), (b), (c) , substitute, into Equation A. 

(K k + K c ) t+At e = t+At Q + t+A V 

The nodal point temperatures in each element can be found by using 
Equations (la), (lb), (lc). 

II) Nonlinear Steady-State Conditions 

The governing equation of nonlinear steady-state heat transfer is 

[ fc K k + t K C + t K r ] A6 (l) = t + At Qt t + At Q c (i-1) 

+ t - At Q r (i - 1 ) - t + At Q k (i - 1 ) 
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where 
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1 ♦ “q 


Vector of concentrated nodal point heat flow input 


k = Thermal conductivity 
h = Convection coefficient 
a = Radiation coefficient 


a 


5.667 x 10 
5.667 x 10 


ef (®r * 0 « )( 6 r * 8 «) 

= Stefan - Boltzman Constant 


c = Emissivity 
f = Viewfactor 
0^ = Surface temperature 
0^ = Temperature of radiated body 


The thermal conductivity, k, and radiation emissivity, £, are temperature 
dependent. The convection coefficients, h, are both time and temperature 
dependent. Because of this, an iterative procedure must be used to solve 
these equations. A Newton-Raphson method is used to iterate to convergence 
the following equation: 


t + At 0 (i) _ t - At 0 (i - 1) + t + At AQ (i) 
Convergence is evaluated by both of the following: 


I e (l) - e (l " | ^ 6. 


a (i) - a (i " 

- “ 1 < * 

.(i - 1) ' = 6 2 


The coefficients k, h, and a are input data. They are set up as a data base 
(or data table) and stored on a file for use during the iterative procedure. 


39 



Ill Transient Analysis 


In transient heat transfer analysis the heat capacity effects must be 
included in the analysis as part of the rate of heat generated. If the Euler 
backward implicit time integration is employed the heat flow equilibrium 
equations used are obtained directly from the equations of steady-state 
conditions; 


0^ (x,y,z) = (x,y ,z) 6 (t) 


then 


t+At _ r jj(«)T { t+At b(m) _ t+At c (m) fi (m) . t+At - ) dy 
m v 

whcr€ ^ ^ no longer includes the rate at which heat is stored within the 

material. Thus, the transient conditions are, in linear analysis, 


c t+At § + (K k + k c ) t+At e = t+At Q + t+At Q e 

and in nonlinear analysis, 

t+At q( 1) t+At g(i) + j-t^k + t R e + t R r^ £ q^ 1 ) 
= t+At Q + t+At Q c(i-l) + t+At^rCi-l) _ t+At Q k(i-l) 

where C, t+At C*' 1 ^ are the heat capacity matrices, 


C= £ / h« T . C W H W . dv W 

m m 
v = 

t+At^(i) p jj(m)T t+At c (m)(i) g(m) dy (m) 

L m y(m) 
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Graded-Composite Structures 


The structure is built up from a series of layers of different materials, 
such that the material properties are discontinuous functions of an appro- 
priate integration through the thickness has to be carried out. For the CSTEM 
finite element code we use a layered approach, wherein a midpoint rule inte- 
gration scheme is adopted for each layer. 

Layers are numbered sequentially, starting at the bottom surface of 
the shell. Each layer contains stress points (or temperature) on its mid- 
surface. * The stress components of the layer are computed at these stress 
points and are assumed to be constant over the thickness of each layer, so 
that the actual stress distribution of the shell is modelled by a piecewise 
constant approximation as shown in Figure B-3. Layers of different thickness 
can be employed, as well as different number of layers per element. 


Number of Layers per Element 
+h/2 +1.0 



Diagram 

Figure B-3. Layered Model and the Corresponding Stress (Temperature) 
Representation . 
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The matrices for the heat transfer equations are obtained by the inte- 
gration of each layer in the element. 

K k = I / B (i)T t K (i) . B (i) dv (i) 

i=l v 

t+At k(j-l) _ y S B ii)T t+At Q (j-l)] dv (i) 

^ i-1 v 

c = h / H (i)T . c (i) . H (i) dv (i) 
i=l ** v — 

where i is number of the layer and j is iteration point number. 

The thermal conductivity matrix [k] is written as follows to handle the 
orthotropic material properties at each layer of the element. 




Figure B-4. Inclination of the Principal Axes 
of Anisotropy, 
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Computer Program - Flow Charts 


The present computer code considers the graded composite structures. It 
will handle the linear and nonlinear, steady-state, and transient heat 
transfer problems by using the "COLSOL" solution technique and the Modified 
Newton-Raphson iteration scheme to solve the graded composite structures. 

The computer program flow chart is as follows: 
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Nonlinear Transient 
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HEAT TRANSFER 
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CSTEM INPUT SHEETS 


I. HEADING AND CONTROL INFORMATION . 

1.1 Title Card 

1.2 Problem Size Data 

II. NODAL COORDIN ATE S 

IV. ELEMENT DEFINITIONS 

IV. 1 Header Line For Each Element Type 
IV. 14 8 Noded Solid 
IV. IS 16 Noded Solid 
IV. 16 20 Noded Solid 

V. ELEMENT SKEW MATERIAL COORDINATE SYSTEMS 

VI. ELEMENT LAYER DEFINITIONS 

VI. SKEW COORDINATE SYSTEM DEFINITIONS 
VIII. FIXED BOUNDARY CONDITIONS 

HEAT TRANSFER INPUT DATA 

************************************************** 

1) . SUBTITLE 

1.2). LOAD CASE INFORMATIONS 

2) . LOAD CASES CONTROL FUG 
2B). TIME STEPS 

2C). TEMPERATURE STEPS 

3) . THERMAL CONDUCTIVITIES 

4) . INTERNAL HEAT GENERATION 

5) . INITIAL CONDITIONS 

6) . CONVECTION BOUNDARY CONDITIONS 
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7) . HEAT FLUX AT THE SURFACE OF THE BODY 

8) . NODAL POINT OF HEAT FLOW 

9) . PRESCRIBED NODAL TEHPERATURES 

10) . SPECIFIC HEAT 

11) . RADIATION BOUNDARY CONDITIONS 


I. HEADING AND CONTROL INFORMATION 

1.1 Title Card 
IDENTIFICATION 

Up to 60 characters for Identification 

1.2 Problem Size Data 

NN HMAT NLC INGEOH IHTR NDATCK 


NN number of nodes In the model 

MMAT nimiber of different materials In the model 
NLC nimiber of load cases 

INGEOH geometry generation flag 

■ 0 geometry read In on Input file 

■ 1 geometry pregenerated 
(contained on 3 pregenerated files: 

node-element, transformation, material layer) 
IHTR HEAT TRANSF ER OPTIONS: 

■ 1 linear steady-state analysis 

■ 2 nonlin ear steady- state analysis 
- 3 linear transient analysis 

■ 4 nonlinear transient analysis 

(English units used If IHTR >0, metric units 
used If IHTR < 0) 

NDATCK -data check flag 0 < data check only 
> 1-run the problem 


II. NODAL COORDINATES 


(enter NN lines Of 1L. only If INGEOH equals zero) 
N X Y Z 


N node number 

X X coordinate of the node ( in. or cm.) 
Y Y coordinate of the node ( In. or cm.) 
Z Z coordinate of the node ( In. or cm.) 
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IV. ELEMENT DEFINITION 


(enter NELTYP of the following) 

IV. 1 Header Line For Each Element Type 

NTYPEL HELEHS INCOHP INTORD 

NTYPEL element type; 8-8 nodes; 16-16 nodes; 

20-20 nodes 

NELEMS number of elements In this group 

(not relevant when using pregenerated geometry) 
INCOHP 0 - Include Incompatible modes (If applicable ) 
1 - do not Include Incompatible modes 
INTORD Integration order ( If applicable 1e. 2 or 3) 


(enter element definition lines IV. 2 - IV. 16 only If no 
generated or pregenerated geometry (INGE0H-O)) 


IV. 14 8 Noded Solid (NTYPEL - 81 

NEL N1 N2 N3 N4 N5 N6 N7 N8 I MAT 

NEL el ement number 

N1--N8 nodes defining the element, labeled as In Flg.l 

IHAT material number (-1 If element removed) 

IV. 15 16 Noded Solid f NTYPEL - 161 

NEL N1 N2 N3 N4 N5 N6 N7 N8 

N9 N10 Nil N12 N13 N14 N15 N16 IHAT 

NEL el ement number 

N1--N16 nodes defining the element, labeled as in Fig. 2 

IHAT material number (-1 If element removed) 

IV. 16 20 Noded Solid (NTYPEL - 201 

NEL N1 N2 N3 N4 N5 N6 N7 N8 N9 N1Q 

Nil N12 N13 N14 N15 N16 N17 N18 N19 N20 IHAT 

NEL el ement number 

Ml— N20 nodes defining the element, labeled as In Fig. 3 

IHAT material number (-1 If element removed) 


V. ELEMENT SKEW MATERIAL COORDINATE SYSTEMS 

(enter JL only If IN6E0N equals zero 
If no skew , Input 0, then go to VI) 

(axes define principal directions for orthotropic materials) 

i 

NESKEW -number of element skew coordinate systems 

defined* Enter NESKEW of the following lines. 


NEL NTRN NEND N1NC 


NEL element number having a local material coordinate 

system __£znri ' 

NTRN transformation number defining this skew 
NEND last element having this local coordinate system 

NINC Increment on elements from NEL to NEND having 

this local coordinate system 
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VI. ELEMEHT LAYER DEFINITION 

(enter VI. 1 only If INGEOM equals zero; IF not composite 
structures Input 0, then go to VII) 

VI. 1 Manual Laver Input 

NELAY 


NELAY number of layer definitions to be read 
(enter NELAY of the following line sets) 


NEL NLAYR LAX NEND NINC 


NEL element for which this layer definition 

applies 

NLAYR number of layers In this element 
LAX local element system axis perpendicular to 

layers (1,2, or 3) 

NEND ending element with this layer system 
NINC Increment In element number 

(enter NLAYR of the following line) 

(enter layers In the order as encountered when moving 
In + direction along LAX) 


I MAT IP^R NTRN ANS 


IMAT material number of this layer 

I PER Integer percentage of element thickness 

occupied by this layer 

NTRN transformati >n number of skew system, 

referencing orientation of material axes 
with respect to global (If NTRN»0, ANG Is 
used) 

ANG angle rotated about element axis LAX, 

referencing orientation of material axes with 
respect to element local (used only If NTRN-0) 

NOTE: Rotations are always positive 

counter-clockwise as viewed from the positive 
end of the axis of rotation. 

If rotating about X, 0 Is referenced to the Y axis. 

If rotating about Y, 0 Is referenced to the Z axis. 

If rotating about Z, 0 Is referenced to the X axis. 



VII. SKEW COORD I MATE SYSTEM (NTRN REFER TO THIS SYSTEM) 
VH.l Reference Nodes 

NSKEWR (IF NO, INPUT 0, THEN 60 TO VII. 2) 

NSKEWR- number of skew coordinate systems defined with 
reference nodes. Enter NSKEWR of the following 
lines. 


NTRN NO N1 NJ 


NTRN transformation number associated with this 
skew system 

NO node on local x axis 

NI node on x axis In +x direction 

NJ node In xy plane In general direction of +y 

axis 


VII. 2 Direc tion Cosines 


NSKEWD (IF NO, INPUT 0 , THEN 60 TO VII. 3) 


NSKEWD- number of skew coordinate systems defined with 
direction cosines. Enter NSKEWD of the 
following lines. 


NTRN All A12 A13 A21 A22 A23 


NTRN transformation number associated with this 
sic6w sy$tM 

All, A12, A13 - direction cosines defining local x axis. 
The local X axis Is defined by ■ All*i + 
A12*i ♦ A13*k 

A21, A22, A23 ■ direction cosines defining local y axis. 
The local y axis Is defined by §y ■ A21*l + 
A22*i + A23*fc. The local 2 axis Is In the 
direction such that It Is a right handed 
coordinate system. 
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VI I. 3 Euler Angles 


NSKEWE (IF NO, INPUT 0, THEN 60 TO X4) 

NSKEWE- number of skew coordinate systems defined with Eul 
angles. Enter NSKEWE of the following lines 

» 

NTRN ALPHA BETA 6AHHA 

NTRN transformation number associated with this skew system 
ALPHA, BETA, GANNA Euler angle successive right 

hand rotations along the 
local x, y, i z axes, beginning 
with the global system 


(enter only If using ICAN routines to generate elastic 
material properties) 

X. 1 ICAN GENERATED ELASTIC PROPERTIES 

NL NNS 


NL number of layers 

NNS number of different material systems 

l Ente* NL of the following cards.) 

INP1 IP1 TU TCU OELN 


INP1 layer number (Identifier for the layer) 

I PI material number (which of the mat'l cards 

applies to this layer) 

TU use temperature (used for linear steady state 

heat transfet analysis) 

TCU cure temperature 



(Enter MMS of the following lino) 

ci c 2 vfp wp C3 £4 55£ m — m 

Cl primary fiber type data bank code 

* 

C2 primary matrix type data bank code 

VFP primary fiber volume ratio 

WP primary void volume ratio 

C3 secondary fiber type data bank code (same 
as Cl for standard composite systems). 

C4 secondary matrix type data bank code (same 
as C2 for standard composite systems). 

VSC secondary composite system volume ratio 
(-0 for standard composite systems) 

VFS secondary fiber volume ratio 

WS secondary void volume ratio 

NOTE: 


Material numbers are assigned consecutively In order 
of Input, l.e. first material system Input Is material 
number 1, second material system Input Is material 
number 2, etc. 
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HEAT TRANSFER INPUT DATA 


(REPEAT NLC TIMES FROM ITEM 1 TO THE END OF LOADING 
CONDITIONS) 



IHTR -1 linear steady-state 

■2 nonlinear steady-state 
■3 linear transient 
■4 nonlinear transient 

NICAN -1 use ICAN to generate elastic properties 
-O no use ICAN 



NOTES: 


NHGEN ( 0- no; 1-yes) Internal heat generation 

NCBC ( 0- no; 1-yes) convection boundary conditions. 

NFBC ( 0- no; 1-yes) heat flux at surface boundary 

conditions. 

NITEM ( 0- N0;l-y.s) Initial condition 

NRBC ( 0- no; 1-yes) radiation conditions 

NODFLW ( 0- no; 1-yes) nodal point heat flow Input. 

ITER The max. Iterations. 

CDELTA Convergence criteria; 

ABS 1 1 [ T( 1 ) -T ( 1 - 1 ) ]/ T(1) || < CDELTA (l.e. 0.01 ). 
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NT I MSP -number of tlM load cases 

TIMSEP(NTIMSP)- time at each time segment, (unit-hour) 
TIMINC(NTIMSP-l) tlie Increment between the time 

steps. 


2c). TEMPERATURE FUNCTIONS (temp, step) 

(only for nonlinear ANALYSIS-2 A 4 ; OTHERS skip this line) 



NOTE: 

NTEMP number of the temperature points for the temp, 

dependent curves 

TEMPSP(NTEMP) Temperature at each temp, point. (degree C or F) 



NOTES: 


MMAT NO. OF MATERIAL TYPES IN THE STRUCTURE. 

MMT MATERIAL NUMBER 

KXX( MMT, NTEMP) | 

KYY( MMT ! NTEMP) ) Thermal conductivity table In the 
KZZ(MMT, NTEMP) | PRINCIPAL AXIS. 

(H/M-C) or (8TU/HR-FT-F) 
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4a). INTERNAL HEAT GEHERATIOH flf NHGEN-O; INPUT 0, GO TO 5) 
NEL. NSQ. NEND. INC 


9 (LAST L I NE ) 


4bl . DATA SET OF HEAT SOUUCE (TINE DEPENDENT). 

NSEIHfl 

HQAT (I.!) HQAT ( 1 « NTIMSP ) 


HQATfNSETHQ.n HQAT (NSETHQ. NT I MSP 1 

NOTES: 

NEL The beginning element # In this group. 

NSQ The data set # of this group. 

NSETHQ The total number of data set In the Internal 
heat generator 

HQAT (NSETHQ, NT I HSP) - Heat generated rate per volume. 

( Nat/ [meter ]**3; or BTU/[hr-ft**3] ). 
at each data set of each time step 
NEND End of the element In this group. 

INC Increment of the element In this group. 


5). INITIAL CONDITION 

(If NITEH .eq. 0, Input 0, go to 6) 

NN. TEMP. NEND. INC 


0 (LAST LINE) 

NOTES: 

NN - The beginning node number # In this group. 

TEHP - Initial nodal temperature, (degree of F; 

or degree of C). 

NEND - End of the node In this group. 

INC - Increment of the node In this group. 


6a). CONVECTION BOUNDARY CONDITION . 

(if NCBC .eq. 0; Input 0 and go to 7). 

NEL. I FACE. NDS.NEHO, INC 

OL a sLiim 

6b). DATA SET OF CONVECTION TENF. AT THE ENVIROHENTS 
NSETCO 

TEHPCOn . U TENPCQTUfUISPl 

TEHPCOTNSETCO 1 1)1 1!!! I !TE]ipCO(NSETCO.NTIHSP) 

NOTES: 

NEL The beginning element # in this group. 

IFACE The face i In this element. (the definition 

see the figure). 

NDS The data set # in this group. 

NEND The end of element In this group 

INC Increment of the element In this group. 

NSETCO The data set of convection temperatures. 

TEHPCO(NSETCO.NTINSP) temperatures at each data set 
(NSETCO) .and each time step (NTIHSP). 

(DEGREE OF C OR F) 
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6c). CONVECTION COEFFICIENT TABLES % 

(Enter NSETCO convection coefficient tables) 

HCOV(l,l) HCOV(l,NTINSP)) 


HCOV(NSETCO, 1) ,HCOV(NSETCO,NTIMSP) 


NOTES: 

NHT Material number 

NTEMP No of temp, points In the coeff. curve 

HCOV (NTEMP, MMAT) Coeff. of convection, 

of each material at each temperature 
point (temp dependent) 

( Wat/[meter**2-degree C] or BTU/[ hr-ft**2-degree F] ). 


7a). HEAT FLUX AT THE SURFACE 

( If NFBC .eq. 0 then Input 0 , go to item 8) 

NEL. I FACE. NSF. NEND. INC 

6(LAST’liNEj* 

7bl.HEAT FLUX TEMPERATURES. 

NSETFl 

TEMPFU1.1) TEMPFUl.NTIMSPl 

• • • i 

T EMPFLf NSETFL. 1) TEMPFUNS ETF l.NTIMSP) 


NOTES: 

NEL The beginning element # In this group. 

IFACE Heat flux at the surface # In the element NEL. 

(the face definition see the figure.) 

NSF The data set # in this group. 

NSETFL Total dat set in this group. 

TEMPFL(NSETFL.MTIMSP) - Heat Flux at the surface, 
at each data set and each time step. 
(Wat/[meter**2]; or BTU/[ hr-ft**2] ). 

NEND The end of the element In this group. 

INC Increment of the element in this group. 


8a). NODAL POINT HEAT FLUX 

( If NODFLH .eg. 0; Input 0, go to Item 9) 

Nl. NSP. NEND. INC 
6 (LAST’ LINE) 

8b). DATA SET OF NODAL TEMPERATURES 
NSETPQ 

PQ(I.l) > «»»»»» u 1 1 1 • « PQf 1 .NTIHSP) 


PQ(NSETPQ.l) ............ PQfNSETPQ. NTIHSP) 

NOTES: ... ■- 

Nl The beginning node # In this group. 

NSP The data set # In this group. - 

NSETPQ Total data set In this group. 

PQ(NSETPQ,NTIMSP) Heat flux nodal point Input, 
at each data set and each time step 
( Hat. or VTU/hr ). 

NEND The end of node # In this group. 

INC Increment of the node # In this group. 
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9a). PRESCRIBED NODAL TEMPERATURES 

(THIS ITEM MUST BE INPUT; IT LIKES STRUCTURAL SUPPORT, 
WITHOUT THIS INPUT IT WILL CREAT A RIDGE-BODY NOTION) 

NN. NSPR. NEND. INC 


OILAST LINE) 


9b). TEMPERATURES AT EACH TIME STEP 

NS ETP B 

PRESTM(l.l) ............... .PRESTMfl.NTIMSPl 


PRESTMINSETPR.l) PRESTM(NSETPR.NTIMSP) 

NOTES: 


NN The beginning node # In the group. 

NSPR The data set # In this group. 

NSETPR TOTAL DAT SET on this group. 

PRESTM (NSETPR , NT I MSP ) - The prescribed nodal 

temperature, (degree F or C ). 
at each data set and each time step. 
NEND The end of the node In this group. 

INC Increment of the node In this group. 


10). SPECIFIC HEAT (only for transient analysis). 

(If steady-stata analysis skip this line go to 11). 


MASSY 1 ) HASS(MHAT) 

HSPCn.l) HSPC(l.NTEHP) 


HSPC(MMAT.l) HSPC(MHAT.NTEMP) 

NOTES: 

MASS(MMAT) Mass density (Kg or Ib.m) at the 

material MMAT 

HSPC(MMAT.NTEMP) THE SPECIFIC HEAT of the material 

MMAT, at each temp, point. 
(KJ/kg-C) or (BTU/lbm-F) 


11a). RADIATION BOUNDARY CONDITIONS 

(If IANALY .eq. 2 or IANALY .eq. 4 Input this line; 
otherwise skip 11, 60 TO ITEM 1 FOR NEXT LOAD CASE) 
(If NRBC .eq. 0 Input 0) 

NEL. IFACE. NSR. NEND. INC 


0(LAST LINE) 


llb>. RADIATION TEMPERATURES. 


NSETRA 

TEHPRAfl.I1 TEHPRAfl.NTIHSPI 

TEMPRAf NSETRA. 11 ....... .TENPRAf NSETRA. NTIHSP1 


NOTES; 

NEL The beginning element # In this group. 

IFACE The face of the element which the radiation 

TEHPRA(NSETRA,NT INSP) - The temperatures at the 

surface nodes. In the data set and each time 
step, (degree of F or C ). 

NSR The data set # In this group. 

NSETRA Total data set In the group. 

NEND The end of the element In this group. 

INC Increment of the element In this group. 


llcl. VIEWFACTOR 


HEi VPU W .NE NPi N IM£ 


0 (LAST LINE) 

NOTES: 

NEL # OF ELEMENT 

P(6) THE VIEW FACTOR AT THE FACE OF ELEMENT 
NEND LAST ELEMENT OF THIS 6R0UP 

NINC INCREMENT OF THE ELEMENT 
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lid). EMISSIVITY TABLES - (enter MHAT times). 

EHCOV(HHT.l). EHC0V(HHT,2) EHCOVfHHT.HTEHP) 


NOTES: 


EHCOY(MNT,NTEMP) The emlsslvity coeff. at each temp, point 
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Linear Heat Transfer - Composite Structures. 
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REFERENCE NODES 1 1 0 1 1 0 0 200 0.010000 

****** THE TINE FUNCTION ****** 

NSKEWR THE TOTAL NUMBERS OF TIHE SEGMENTAL POINTS - 6 




CONVECTION NODAL TEMPERATURES 
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MATERIAL 1 THERMAL CONDUCTIVITY TABLE 
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LINEAR TRANSIENT HEAT TRANSFER 



LINEAR STEADY-STATE HEAT TRANSFER 
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ELEMENT TYPE 20 DIRECTION COSINES 

20 NODEO ISOPARAMETRIC BRICKS 

HSKEWO 




REFERENCE NODES 1 1 0 1 1 0 0 200 0.010000 

****** THE TIME FUNCTION ****** 

NSKEWR V ' THE TOTAL NUMBERS OF TIME SEGMENTAL POINTS - G 


CONVECTION NODAL TEMPERATURES 



PRESCRIBED NODAL TEHPERATURES 
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